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We obtain a complete set of free-field solutions of the Dirac equation in a (longitudinal) boost- 
invariant geometry with azimuthal symmetry and use these solutions to perform the canonical 
quantization of a free Dirac field of mass M. This coordinate system which uses the 1+1 dimensional 
fluid rapidity r\ = l/21n[(t — z)/(t + z)] and the fluid proper time r = (t 2 — z 2 ) 1 ^ 2 is relevant for 
understanding particle production of quarks and antiquarks following an ultrarelativistic collision of 
heavy ions, as it incorporates the (approximate) longitudinal "boost invariance" of the distribution 
of outgoing particles. We compare two approaches to solving the Dirac equation in curvilinear 
coordinates, one directly using Vierbeins, and one using a "diagonal" Vierbein representation. 

PACS numbers: 25.75.-q, 04.60.Ds 
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I. INTRODUCTION 

The kinematics of heavy-ion collisions at rclativistic 
colliders is such that the distribution of secondary par- 
ticles obeys an approximate longitudinal boost invari- 
ance. This approximate invariance was the basis for 
many hydrodynamical simulations of particle produc- 
tion following the collision of two ultrarelativistic heavy 
ions. The use of this symmetry for heavy-ion colli- 
sions was advocated by Bjorken [l[ in a seminal paper 
which was based on earlier work using Landau's hydro- 
dynamic model [H by Cooper, Frye and Schonberg Q 
and Hwa 0]. In a hydrodynamic description of the time 
evolution of the plasma produced following a collision, 
this invariance leads to the surfaces of constant energy 
density being functions primarily of the kinematic vari- 
able r = (t 2 



j 2 ) 1 / 2 , where the z axis is the axis of 



cylindrical symmetry. For head-on collisions of identical 
nuclei, the kinematics also has cylindrical symmetry. 

In trying to understand the particle production of 
hadrons in this scenario from first principles, one ap- 
proach is to assume that a significant fraction of the ini- 
tial energy that is eventually transformed into the final 
hadrons is in a classical gluonic field configuration which 
produces quarks and quantum gluons by the Schwingcr 
tunneling mechanism [5[ . In this approximation the par- 
ticles are produced by nonperturbative tunneling. How- 
ever, we ignore production due to further quantum cor- 
rections due to rcscattcring of the quarks and gluons that 
are "sparked" from the vacuum. 

In this boost-invariant situation, all physical quantities 
such as energy densities are just functions of the proper 
time t. In particular the chromoclectric field density E a , 
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a = 1, 2, ... 8 is also a function of r. Choosing the chro- 
moelectric field in the longitudinal direction, then we can 
choose a gauge where the vector potential is also a func- 
tion only of u = In r and is defined via the relationship 



El(u) = 



dAj{u) 
du 



(1.1) 
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In this picture the "classical gluonic field" would ap- 
proximately be a function of r and one would obtain the 
initial quantum plasma by studying the coupled problem 
of the produced quarks and gluons interacting with the 
background classical gluonic field. Such a picture was 
advocated by Cooper et al. [f| for the production of par- 
ticles in the central rapidity region. Once such a system 
of equations is solved, one then has to construct the phys- 
ical quantities of QCD which depend on the two Casimir 
invariants of color SU(3). It is only these quantities which 
are gauge invariant. 

In the hydrodynamic scaling regime where the fluid ve- 
locity v = z/t, the fluid rapidity 77 = l/2ln[(t-z)/(t+z)] 
can be identified with the particle rapidity: l/21n[(l — 
f)/(l + v)]. In Ref. [|| numerical simulations were per- 
formed in 1+1 dimensions, with an a Abclian external 
field. Although the formalism for doing the 3+1 dimen- 
sional calculation was presented, simulations were not 
done in that paper and the cylindrical symmetry of the 
head-on collision problem was not exploited. By exploit- 
ing the cylindrical symmetry, the computational com- 
plexity of the back reaction problem is reduced from N 3 
modes to describe the quantum field of the quarks to N 2 
modes which is quite a savings. 

Recently the problem of finding the transverse mo- 
mentum distribution of quarks and gluons produces by 
a strong constant chromoelectric field by the Schwingcr 
mechanism [j| has been solved Q ■ This calculation found 
the interesting physical result that the transverse distri- 
butions of produced particles depends not only on an in- 
variant related to the energy density of the plasma, but 
also an invariant which describes the direction in color 



2 



space that the initial chromoelectric field was pointing. 
This might have important experimental consequences. 
In this paper we will discuss the canonical quantization 
of the free Dirac field in the coordinate system which 
utilizes the curvilinear coordinates r and the "fluid" lon- 
gitudinal rapidity 77 = | ln[(i — z)/(t + z)] as well as the 
cylindrical coordinates p and <f> that correspond to the 
directions transverse to the collision axis. This formu- 
lation will simplify identifying the adiabatic solutions to 
the back reaction problem and will simplify the computa- 
tional effort needed to determine transverse momentum 
distributions. The canonical quantization of the fermion 
field presented here will also allow us to precisely define 
the initial state of the quarks before the chromoelectric 
field is "turned on" , so that the initial values of the time 
evolution problem and back reaction can be specified. 

The problem of quantizing the Dirac field in an arbi- 
trary coordinate system was first discussed in the classic 
paper of Brill and Wheeler 0] who were concerned with 
understanding the motion of electrons and neutrinos in 
a background gravitational field. This problem was also 
later considered by Parker Q and collaborators who were 
interested in the problem of particle creation in expand- 
ing universes. More recently, Viallalba and coworkers 
have concerned themselves with the general problem of 
separation of variables [10| in external vector fields in or- 
der to find exact solutions [ll| of the Dirac equation. The 
problem of quantizing the Dirac equation on the hyper- 



boloids 



■y 



-z 2 ) 1 / 2 was also studied in g reat 



detail in the I970's by Sommerfield [12J, diSessa [13J and 
others Furthermore, the problem of quantizing the 
Dirac equation in cylindrical coordinates was studied by 
Balantekin and DeWeerd [l5| also with interest in ultra- 
relativistic heavy-ion collisions, but without exploiting 
the boost invariance aspect of this problem. 

The intent of this paper is to fill the apparent gap in 
the previous work on this subject by studying the Dirac 
equation in a situation where the 1+1 dimensional fluid 
proper time r = (t 2 — z 2 ) 1 / 2 is important as well as the 
cylindrical symmetry. In doing so we will also compare 
and relate two different ways of conside ring this problem, 
first directly using Vierbeins (tetrads) [16[ to transform 
into the curvilinear coordinate system, and the second 
making a similarity transform into a frame where the 
Vierbeins are diagonal. 

The paper is organized as follows: In Sec. |TT] we will 
review the general mathematical framework of the sub- 
sequent discussion and outline our approach. In Sec. IIIII 
we will discuss the 1+1 realization of this problem using 
boost-invariant coordinates. In Sec. IIVI we will discuss 
the case of a 3+1 dimensional Dirac equation in cylindri- 
cal coordinates. Finally, in Sec. [V] we will combine the 
results of the previous two sections in order to discuss 
the solution of the Dirac equation in longitudinal boost- 
invariant coordinates with azimuthal symmetry. We con- 
clude in Sec. IVII 



II. GENERAL FRAMEWORK 

The free-field Dirac equation in Cartesian coordinates 
(£") is given by 



( 7 Q d a + M) = . 
Throughout this paper we use the metric 
= diag(-l, 1,1,1), 



(2.1) 



(2.2) 



and the Dirac gamma matrices in the chiral representa- 
tion, in 3+1 dimensions, given as 



7 



r 
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Ok 


-1 




0~k 






(2.3) 



where a k arc the usual Pauli matrices, and 1 is the unit 
(2 x 2) matrix. The constant Dirac matrices satisfy the 
relation 



{7",/} = 277 



a/3 



(2.4) 



Next, we will consider the change of variables 
(£ Q — > x a ). We introduce the Vierbeins, V a a , and their 
inverses, V" as 



IT - £ 

n d dx a d 



Q£a Q£* Q x a 



V a a dx a 
V a d 



We have 



'Jap 



v a Q v Q b = s a 



(2.5) 
(2.6) 

(2.7) 



Under the change of variables the line element transforms 
as 



ds 2 = Vaf3 dCd^ = g ab dx a dx b 



(2. 



In the following we will use the notation y/—g = det[g a f,]. 
One also has the connection between the Cartesian met- 
ric and the metric in the curvilinear coordinates 



gab = V Q a V\ r, af} 



(2.9) 



Finally, in the new system of coordinates the free-field 
Dirac equation (|2.ip becomes 



(Td a + M) 4, 







(2.10) 



To simplify the Dirac equation in an arbitrary curvi- 
linear coordinate system one often performs a similar- 
ity transformation to diagonalize the Vierbein matrices. 
This then produces a representation which we will con- 
note the "diagonal tetrad" representation. To obtain the 
diagonal tetrad representation, one introduces a similar- 
ity transformation, S, such that 



S- 1 7 a S 



r 



(2-11) 
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where the gamma matrices 7° are proportional to the 
constant gamma matrices, 7 . By construction, the two 
set of gamma matrices, 7° and 7 , are defined such that 



7 



V a Q 7 Q 



(2.12) 



where V a a are the Vierbeins in the diagonal tetrad rep- 
resentation. The two sets of gamma matrices, 7° and 7° 
satisfy the same algebra 



(7 , 7 ) = {7 , 7 ) = 25 



(2.13) 



Hence, the similarity transformation S which connects 
these two irreducible sets of matrices, is unique [l7l ]. 

In the diagonal tetrad representation the Dirac equa- 
tion becomes 



(7 a V Q + M) ^ [d] = 



(2.14) 



where we have introduced the diagonal Dirac wave func- 
tion 



-0[d] = S 1 ip , 
and the covariant derivatives 

V Q = S- 1 d a S = d a - T a 

Here, the spinor connection T a is given by 



1 



a) > 



(2.15) 
(2.16) 

(2.17) 



and we have introduced the notation E a/3 = j[y a , 7^], 
defined in terms of the constant Dirac matrices. 

The Christoffcl symbols, r£ , are obtained from the 
Lagrange equations of motion 



d_ / d£_ \ _ dC_ 



At \dx a 



dx a 



= 0, 



(2.18) 



and the conservation law can also be written as 

J A d+1 x M (</=s V") = 0. (2.23) 

or 



dE^ y/=g V M = . 



(2.24) 



Here E M is the rf-dimensional hyper-surface perpendicular 
to the \i direction. Choosing \i = then the conserved 
charge is (with V M = j M = V'7 M ip) 



J xq— const. 

Thus the natural Dirac inner product is : 



(-01, "02 



A a x 



-9 (^1 7° "02 



(2.25) 



(2.26) 



with ip = ip^j . The canonical quantization of the Dirac 
field assumes a mode expansion: 



ip{x) 



A d k 



^l + \x)a k + ^\x)b 



with 



{a k , at,} = {6 k , bl} = 5(k'-k) 



(2.27) 



(2.28) 



Here (x) are positive- and negative-energy solu- 
tions [20J of the homogeneous Dirac equation and satisfy 
the orthonormality relations 



(4^ i 



= £(k' - k) , 
. 



(2.29) 
(2.30) 



In the following, we will study the canonical quantiza- 
tion of the free Dirac field in the Cartesian and diagonal 
tetrad representations, in the context of several system 
of coordinates characteristic to the physics of heavy-ion 
collisions. 



x a +T a hc x b x c 



. 



(2.19) 



corresponding to the Lagrangian of a free particle in the 
x a coordinates llal. i.e. 



C = g a b i a x b . 



(2.20) 



The inner product is defined by considering the covari- 
ant conservation law for a current. In an arbitrary num- 
ber (d) of spatial dimensions, we have the conservation 
law fijl 



-9 n 



0, 



(2.21) 



provided that V M vanishes at infinity. Here the covariant 
divergence is defined as 



V M = 



1 



(2.22) 



III. BOOST-INVARIANT COORDINATES IN 
1+1 DIMENSIONS 

For illustrative purposes, we consider first the 1+1 di- 
mensional case studied previously by diSessa [l3[. Here, 
the Dirac equation is 



(7°<9o + 7 3 <9 3 + M) ip = 0. 



(3.1) 



In 1+1 dimensions, the chiral representation of the 
Dirac gamma matrices for our choice of metric, r) al3 = 
diag(— 1, 1), is (compare with Eg. 12.3ft 



7 



1 

-1 



7 



1 

1 



Consider the change of variables (£ Q — * x a ) 
t = t cosh r\ , z = r sinh rj , 
with the metric g a b = diag(— 1, t 2 ), and 



(3.2) 
(3.3) 
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A. Cartesian tetrad representation 

For our change of variables, the Vierbeins are calcu- 
lated as 



cosh rj t sinhr/ 
sinhry r cosh 77 



v a 

whereas the inverse Vierbeins are 
V" = 



cosh rj — sinh rj 
— sinh 77 — cosh 77 

T T 



We can now write the Dirac equation as 
(T d a + M) V - , 
where the Dirac matrices are given by 



7 = 7 cosh 77 — 7 3 sinh 77 
1 

T 



7 3 = — (—7° sinh rj + 7 s cosh rj) 
t x ' 

Equation (|3.4|) can be explicitly written as 

= 0, 



(3.4) 

(3.5) 
(3.6) 

(3.7) 

where we have introduced the raising and lowering oper- 
ators 

T_ = d t + d z = e - ' (d T + - d„) , (3.8) 
T+ = d t - d z = e" (d T - i a„) , (3.9) 
which obey the properties 
71 



M 71 " 






-r+ m 




^2 _ 



3 ife "" J ifcl) (Mr)j = +Me^- 1 ^J ikri . 1 (Mr), 

(3.10) 



e ife ""J ife „(Mr) 



M e (ik n+ i) v j ifet;+1 (Mr) , 

(3.11) 



where Jifc (Mr) are Bessel functions of complex order 
(and real argument) plj . Throughout this paper we use 
the definitions and properties of the Bessel functions as 
presented in Ref. |22j . 

Therefore, we find the amplitudes ipi and ^2 obey the 
equations 



(71 r+ + m 2 ) ^1 = 0, 
(r+r_ + m 2 ) v 2 = . 



(3.12) 
(3.13) 



Alternatively we could have "squared" the Dirac equa- 
tion (|3.4p in the usual textbook fashion. This would have 
doubled the number of solutions and we would have kept 
only one set such as done in [23|. 

The solutions of the above equations have the general 
form 



V(r,v7)= / dk v i> kr) {T,r]). (3.14) 
It is necessary [24| to seek a solution ip k of the form 



e-^ J ikri _ k {Mr) 
-e' /2 4 + i(«r) 



(3.15) 



Requiring the usual positivc/ncgative-energy boundary 
conditions 1201, we obtain 



e-^H™ .(Mr) 



, (3.16) 



and 



e -n/2 ff (i) i t Mr s 

1 V 2 

- e"/ 2 if.W , t (Mr) 



(3.17) 



Then, the canonical quantization of the Dirac field is ob- 
tained as 

VKr, ?? ) = J dk v [^{r^a^+^iT,^]. (3.18) 

The normalization constants N k are obtained by re- 
quiring the fields V'i ± ' > to satisfy the orthonormality re- 
lations 



(3.19) 
(3.20) 



where the inner product (|2.26[) gives 
^k' v ) = - t dr/ ip kv 7 ip ki) 



(3.21) 



r / drj 



(+)*./.(+) 



with ip = ip' 7 . We work out each term separately: 
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r/,\+) e " = N^*N { + ] i e'^-Mi H ( }\ , t , (3.22) 

4^ e"" = - < ) *^ ) i e^"^ e"^ 2T« +i , (3.23) 

4t 4:1 eV = H) ei(fc;, " fc " ) " ' (3 - 24) 



where we have used the properties of the Hankel functions 
of real argument, x [2^| 

H^{x) = H^ l \x), (3.26) 
H { l?\x) = e^ + -^ v H^ 2 \x) . (3.27) 
We will also use the Wronskian identity [22j : 

H^ ) (x)H^_ 1 (x)-H^ 1 (x)H^(x) = -ii . (3.28) 



such that 



with 



s- 1 r s T = r , 



r-,0 ..0 ^.3 - 1 - ..3 

7=7, 7 = - 7 

T 



(3.35) 



(3.36) 



The diagonal (rotating) Vierbein representation is iden- 
tified as 



We note that by performing the integral in (|3.2ip we 
obtain 



dr] e i{k '- k)ri = 2n S(k' -k) . (3.29) 



Then, we find the normalization constants must satisfy 
the conditions 



1 = M ±} *Ni ±} o Trf ' 



M 



(3.30) 



We pick the normalization constants as: 



^? = \/^^/ 2 . (3.31) 



Therefore, we find that the positive/negative energy so- 
lutions of (|3.4p are 



_ e nk, 1 /2 e ik vV 



e-'/ 2 ^_,(Mr) 



(3.32) 



and 



e - 7rfc,/2 e ifc r) ?7 



^ 2 

-e^tfW (Mr) 

lK.r)T 2 



(3.33) 



B. Diagonal tetrad representation 



We begin with Eq. (|3.4[) and introduce the similarity 
transformation 

S T = cosh | - 7° 7 3 sinh | = exp(-| 7 7 3 ), (3.34) 



V Q Q =diag(l,r), 
V a a = diag(l, -) . 



(3.37) 
(3.38) 



The Dirac equation in the diagonal tetrad representation 
is 



(7 a V Q + M) i> [d] = 0, 



(3.39) 



where we have introduced the diagonal Dirac wavefunc- 
tion 



-0[d] = s T 1 i> . 



(3.40) 



The identification process described in Sec. HI] leads to 
the nonzero Christoffcl symbols 



r T = t 

vv ' 
P?7 _ pn _ _ 



the spinor connections 



1 



r„=2 7°r. r T = o, 



and the covariant derivatives 



1 



0.,3 



V r = d T , V„ = d n - - 7 U 7 

Finally, the Dirac equation becomes 

7° (d T + i-) + 7 3 \ d v + M] V[d] = 
and has a solution of the form 



(3.41) 
(3.42) 

(3.43) 

(3.44) 

(3.45) 

(3.46) 
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We find positive- and negative-energy solutions obeying 
the relations (|3.19p . as 



Yk v [d] y g 



and 



^(-) = *k v /2 ik^r, 

i fc,,[d] 



1& V 2 



V 2 



(3.47) 



(3.48) 



Finally, the Dirac fields in the Cartesian and diagonal 
tetrad representations are related via the transforma- 
tion (|3.40p . Indeed, with our choice of gamma matrices, 
we obtain 



7°7 3 



1 

-1 



(3.49) 



which gives the explicit form of the similarity transfor- 
mation ([3734]) 



S T = 



e-"/ 2 
e"/ 2 



(3.50) 



We immediately verify that the Dirac solutions (|3.32|) 
and (|3.33[) are obtained from (|3.47[) and (|3.48[) as 



,/.(±) _ c ,/.(±) 



(3.51) 



IV. CYLINDRICAL COORDINATES IN 3+1 
DIMENSIONS 



The canonical quantization of the Dirac equation in 
cylindrical coordinates has been studied in detail by Bal- 
antekin and DeWeerd [l5| using the eigenstates of the 
transverse helicity operators to obtain a complete set of 
solutions. Here, for uniformity of presentation we will 
use the tetrad formalism to obtain a related decompo- 
sition in terms of Bessel functions of 1/2 integer order. 
Consider the change of variables (£ Q — > x a ) : 



(4.1) 



t = t , z = z , 

x = p cos 6 , y = p sin 6 . 

with the metric g a b = diag(— 1, 1, p 2 , 1), and yj—g 
A. Cartesian tetrad representation 



The Vierbeins in /iied-tetrad representation are calcu- 
lated as 



V Q = 



10 

cos 9 —p sva.9 

sin 9 p cos 9 

1 



whereas the inverse Vierbeins are obtained as 





I 

cos6> sin^ 

-- sin6» - cos6» 

P P 

1 



The Dirac equation becomes 

(Td a + M) i> = 0, 
with the gamma matrices, 7° = V a a 7 a , defined as 



7° =7° 



7 3 



I 1 

I 2 



7 1 cos 6* + 7 2 sin6> , 



1 

P 



-7 1 sin 6* + 7 2 cos( 



(4.2) 

(4.3) 
(4.4) 

(4.5) 



Following Villalba [lfj , we write Eq. (|4~2j) as 
(Ki + /C 2 )$ = 0, $ = i 7 1 7 2 V 
where JCi and /C2 are the first-order differential operators 

Ki(t,z)= (7°a t +7 3 a z )i7V, 

£ 2 (p,6) = (7 1 a p + 7 2 9 e + M)i 7 1 7 2 , 

and 



(4.6) 



(4.7) 
(4.8) 



17V 



ff 3 
<r 3 



(4.9) 



Since the operators K\ and IC2 commute [231, ^a]=0) 
we can write the eigenvalue equations 

/Ci$ = -/C 2 $ = -A$ . (4.10) 

We will find the solution of Eq. ()4.6|) has the form 

(4.11) 



®{t,p,0,z) = J k x dfc ± <P k± (t, p, 9, z) 



where A 2 - M 2 = k\. 

a. JC2 sector. The structure of the amplitudes <fr k± 
is established by solving the eigenvector problem for the 
operator K.2- We have 



(£a - A) $ k ±(p,0) = 



or 



(i 7 2 d x - i 7 1 d y + i 7V M - A) $ fcj 
which gives 



(4.12) 
(4.13) 



r ~ A + 








V- - 




<$>k ± 


■•Hp,0)' 







-A_ 


-v+ 







($>k ± 


■' 2 (p,o) 


= 0, 





V- 


-A+ 









' 3 (P:0) 










-A_ 






' A (P,0) 





(4.14) 
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where we have introduced the notations X± — A ± M. 
The raising and lowering operators are defined as 



V- = d x — id y = e~ 
V+ = d x + id y = e) 
with the properties 
V- 



d p dgj , 

dp + - dg^j , 



(4.15) 
(4.16) 



V 



e ime J m (k±p) 
e ime J m (k±p) 



(4.17) 



kx \e i(m+1)9 J m+1 (k ± p) 



(4.18) 

The system of equations (|4.14[) decouple, and we obtain 

[P+V- + X+X-] $ k±m {p,8) = 0, (4.19) 
[V-V+ + X+X-] <5> k± '' 3[1] {p,9) = 0. (4.20) 

The spinor Q k± (p, 8) has the form 

oo 

^{p,e)= Y, <Hp,o), (4.2i) 



m— — oo 



with 



A k - e^/\J m _ k (kxp) 
B k - e is / 2 J m+ ,(k ± p) 

C k - e- ie / 2 J m _i(fexp) 
D k - e«V 2 J m+h (k ±P ) 



(4.22) 



and obtain the constraints 



D k± = — A k± = — A kj 

k± A_ 



B 



-L c k± = — C kl 



A 



(4.23) 
(4.24) 



We also find 



A+A_ = A 2 - M 2 = k\ , 



(4.25) 



which gives 

A 2 = J]_ = M 2 + k\ , or A = ± uj± . 

Corresponding to the two eigenvalues, A = ± uj±, we 
two linearly-independent eigenvectors 

A k - er*l 2 J m _ h {kxp) 

^C k - e^ 2 J m+i (k ±P ) 
C k - e-^J m _ l2 (k ± p) 
^A k - e iS / 2 J m+i (k^p) 



(4.26) 
have 



^ s (p,e)=e i 



(4.27) 



with A = swi, and s = ±1. 



b. JCi sector. The amplitudes A k± and C k± are de- 
termined by solving the eigenvalue equation for the op- 
erator JCi, i.e. 



(id + A) ^ s (t,z) = 0, 



(4.28) 



i 1 ° 1 1 1 2 d t +irVl'd z +X 



3_,1_,2 



*nt'' 8 (t,z) = 



(4.29) 



Once again using the representation (|2.3p of the Dirac 
gamma matrices, we obtain the anticommuting products 
of gamma matrices: 



i 7 ° 7 1 7 2 



-<7 3 

a 3 



i 7 3 7 1 7 2 



-1 
-1 



(4.30) 



Since dt and d z commute with the above equation, we 
can also Fourier transform the Dirac field in t and z. We 
introduce the field as 



&±< a = I dE e iEt I dk 



-ik z z ^k±;s 



which gives 

~ -03 
0-3 



-1 
-1 



(US) 

Eq. (|4.32|) is equivalent to 

A i(k z - E) 
i(k z +E) X 



k z ;E ' 



c, 



k z ,E 
k ± 
k z ,E 



(4.31) 

= 0. 
(4.32) 

(4.33) 



For a nontrivial solution, we require the determinant to 
be zero, which gives the usual dispersion relations 



E = ± y^ 2 + k 2 = ± \/M 2 + k 2 
We obtain 



r k ± lA A k ± 

k - E ~ k z -E k - E 



<k z + E) k± 
A Ak - E 



(4.34) 



(4.35) 



Hence, the solution to the free-field Dirac equation can 
be written as 

-,e (*. pAz) = i iW , E (*. pAz) (4-36) 

-e-W* J m _ h (k ± p) 
*e«/ 2 J m+ i(k ± p) 



N kZE e 



X-M k f -E ° "m-H 

-^e-^J^fc^) 
x^e ie / 2 J m+ i(fc ±/9 ) 



The canonical quantization of the Dirac field in cylin- 
drical coordinates is obtained as 



ip{t,p,0,z) = 



fkxdkx I dk z jr ]T 

m=-oo s 



(4.37) 







(+) fej_;s 



m,k z 



+ tp. 



tn,k z " rn,k 



k±;s (-) , k ± ;s f 
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where ip m j, y ' correspond to the positive and nega- 
tive energy solutions, respectively. The normalization 
constants N^'f^g, arc obtained by requiring the fields 

ip-Tl to satisfy the orthonormality relations 

(j\k ± M±) _ 6 J£± - k ^ s( jJ _ k)S , Sl 

(4.38) 

«tf +) .^X (_) )=0. (4-39) 
where we have introduced the inner product 

= -/ pd P fdejaz ^7°^£ 

= I pdp I dd I dz ^ ^£ ■ (4 - 40) 

Using the identities 



d6 e 



i(m —m) 



2,7T dm'm — 

2tt 5{k' z - k z ) = I dz e«~ k ^ z 
S(k' x -k±) 



(4.41) 
(4.42) 

pdp J v (k' ± p) J v (k±p) , (4.43) 



we obtain the normalization constant 



M k ± ;s_±_ E - k z j X-M 
k » E - 2n V 2E V 2A 



(4.44) 



for either positive or negative values of E and A. This 
completes our derivation. 



B. Diagonal tetrad representation 

In cylindrical coordinates, we introduce the similarity 
transformation 



S p = cos- - 7 1 7 2 sin- = exp(-- 7 VJ , (4-45) 



such that 



where 



s; 1 T s p = r , 



(4.46) 



and the Dirac equation can be written as 

( 7 a V Q + M) V> [d ] = 0, (4.50) 

where we have introduced the diagonal Dirac wave func- 
tion 



i'ld] = Sp 1 ip. 



(4.51) 



In the cylindrical system of coordinates, the covariant 
derivatives 



V a — S d a Sp — d a T a 



(4.52) 



Here, the spinor connection T a , 



0. 



(4.53) 



are obtained using the (nonzero) Christoffel symbols 



1 p8 ~ 1 Bp ~ - ■ 

The covariant derivatives are obtained as 



(4.54) 
(4.55) 

(4.56) 



With the above definitions, the Dirac equation (|4.50p can 
be written as 



1 



^[d] = • 
(4.57) 



Following the same approach to the separation of vari- 
ables as in the case of the Cartesian tetrad representa- 
tion, with the difference that now the operator /C 2 is 



17V 



(4.58) 



We find that the solution of Eq. (|4.57[) has the form 



1pE[d] 



J k± dk ± J dk z jr t 



m,kz\E[d] 



, (4.59) 



m— — 00 s 



7° -7°, f=7 1 , 7 2 = - 7 2 , 7 3 = 7 3 - (4.47) 

P 

The diagonal Vierbcin representation, 7" = V a Q 7 Q , is 
introduced as 



V a a =diag(l,l,p,l), 
V a a =diag(l,l,i 1), 



(4.48) 
(4.49) 



where the function V'm^-'l: is given by 



(t,p,6,z) (4.60) 
-J m _i(k±p) 

T^TI kl-E ^m+i(^J-P) 

-A J m -±(k±p) 
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where N^'^ is given by Eq. (|4.44[) . The Dirac fields in 
the Cartesian and diagonal tetrad representations (see 
Eqs. (|4.36[) and (|4.60p , respectively), respectively are 
related via the similarity transformation introduced in 
Eq. (|4.51|) . As advertised, we obtain 





- e -ifl/a 











Sp — 





e i0/2 














e -w/2 
















e W/2 



(4.61) 



which leads to 

itl,E = S P ifcl, m ■ (4-62) 



• diagonal Vicrbcins: 

V a a =diag(l,l,p,r), (5.6) 
V\ =diag(l,l,i i), (5.7) 

• gamma matrices in the rotating-tetrad representa- 
tion: 7° = S^ 1 7° S, where 

7° = 7°, 7 2 = -7 2 , f = -7 3 - (5.8) 

P T 

• diagonal Dirac field: 

tf[d] = s- 1 ^ = S- 1 s; 1 i> . (5.9) 



V. BOOST-INVARIANT COORDINATES WITH 
AZIMUTHAL SYMMETRY 



Consider the change of variables (£ Q x a ) 

t = t cosh r\ , z = t sinh rj , 
x = p cos 9 , y = /j sin , 



(5.1) 



which gives the metric g a b = diag(— 1, 1, p , r ), and 

Since the (x, y) and (t, z) subspaces do not interact, 
the transformations we need are indeed very simple. The 
key ingredients are: 

• Cartesian Vicrbeins: 



V" = 



and 



cosh 77 Or sinh 77 

cos 8 —p sin 9 

sin 9 p cos 9 

sinh 7/ Or cosh 77 



cosh 77 








1 





cos^ 

1 

■— sm( 
P 




sin( 



1 



sinh 77 



— sinh 77 






— cosh 77 

T 



gamma matrices in the fixed-tetrad representation: 



7 = 7 cosh 77 — 7 3 sinh 77 , 

7 3 = — (—7° sinh 77 + 7 3 cosh 77) 

T 

7 1 = 7 1 cos 9 + j 2 sin 9 , 

7 2 = - (-7 1 sin6» + 7 2 cos6») . 

P 



(5.2) 
(5.3) 
(5.4) 
(5.5) 



• similarity transformation: S = S T S p . Note that 
the individual similarity transformations commute 
with each other. 



• covariant derivatives: 

V r = d T , V p 

1 . _ 

V 



do 



(5.10) 



A. Cartesian tetrad representation 

The Dirac equation 

{T d a + M) f/> = , 



(5.11) 



is solved using the separation of variables procedure out- 
lined before, see Eqs. ([46]) . ([47]) and (|4~8j) . The mod- 
ified Dirac field <&(r, p, 9, 77) has the form (|4. 1 1[) . and 
we solve the eigenvalue problems defined in Eq. (14.101) . 
Again, since there is no interaction between the (x, y) 
and (t, z) subspaces, it follows that the (p, 9) structure of 
<&(t, p, 9, 77) remains the same as in (|4.27[) . and we only 
need to concern ourselves only with solving the eigen- 
value problem for the operator JC\ in the new set of coor- 
dinates, (t, 77). Explicitly, we seek the amplitudes (r) 



and C^(t) required by ((427)) . 
Consider the eigenvalue problem 

(£1 + A) $ fe ^(r) = 



(5.12) 



Using the raising and lowering operators defined in 
and (|3.8[) . respectively, the above equation reads 



4t(-r)e- W/2 J m -i(k±p) 



- A 





-71 








A 





T+ 


T + 





A 





_ 


-71 





A 



T)e*/ 2 J m+h {k^p) 
(r)e- ie / 2 J m _.(^p) 
At(r)e^J m+ ,(k^p) 



= 0, 



which is equivalent to 



A 



-71 

A 



<Mr) 



(5.13) 



(5.14) 
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The solutions of the above equations are Bessel functions 
of complex order (and real argument) [2lj . So, similarly 
to Eq. (|3.15|) . the Dirac field in the Cartesian tetrad rep- 
resentation becomes 



T^m,k n ^ ' ' m.k r 



(5.15) 

^^- i '')/ 2 J ifcIJ+ ,(Ar)J m+ ,(fc ± p) 
-e-^)/*J.^ +k {\T)J m _ h {k±p) 
J^e^)/* J ikn _ h {\r)J mH {k±p) 

Requiring the usual positive/negative-energy boundary 
conditions [201 , we obtain 



,k ± ;s(+) _ N k±;s i(m6+k n ri) 



(+) 

- e -i(o-iny2 H W_ i{ XT)J m _ h (k ± p) 
-e-^ +i ^Hfl + l(\r)J m ^(k ± p) 

^ ei(0+i „)/ 2jff (2) ; (Ar)Jm+ (k±p) 



(5.16) 



and 



i\jk±;s i(m9+k v r)) 
iV (") 6 



(5.17) 



-e- 1(fl - i " )/2 ^_ 4 (Ar)J m _i(fcxp) 
^e^-^/ 2 H^ + l(Xr)J m+i (k x p) 
-e-^+^H^J h {\r)J m _ k {k x p) 
^ e i(e^)/2 H W_l(XT)J m+h (k x p) 



The canonical quantization of the Dirac field is ob- 
tained as 



/p oo 
k± dk±_ dk n X! 



(5.18) 



m— — oo s 



k± ;s (+) k± :s 



fl, 



k±;s (-) , k±:s f 



k -s (±) 

where tp m k y correspond to the positive and negative 
energy solutions, respectively. The normalization con- 
stants N^ s are obtained by requiring the fields ip^ tkri 
to satisfy the orthonormality relations 

(^■<°(±) ./*i;«'(±h - s ( k '± Z k ^ S (k' -k)S , s, 

(5.19) 

(€x (+) '^;? _) )=°- ( 5 - 2 °) 

where we introduce the inner product 



{^i>%) = -r jp&pj&ejdz f^; K . 

(5.21) 

We find the normalization constants 



k±]S 1 M(X-M) ±nkv/2 
iV (±) _ 2vr V 16A 6 



(5.22) 



B. Diagonal tetrad representation 

In the diagonal tetrad representation the Dirac equa- 
tion takes the form 



+ 7 - de + 7 - d v + M 

p T 

and has a solution of the form 

M = j k ±- dk ±- j dk 1 J2 



(5.23) 

V'ld] = o , 



m— — oo s 



with 



T m.k. 



i 7 1 7 2 $ fcj -jf = 7V fei;s e i (™e+fc>,')) 

' ' YTL *kyt 



- J ifc,-i(Ar) J m _i(fc±p) 
J ifci)+ i (At) J m+ i(k±p) 
-J ifelj+ i(Ar) J m _i(fc_Lp) 
(At) J m+ i(fc_Lp) 



(5.24) 



(5.25) 



Then, the positive- and negative-energy solutions corre- 
spond to 



,k ± ;s(+) _ N k±;s A{ m e+k v v) 
^m,k v [d] (+) 



(5.26) 



f k ] _AXT)j m _ h {k^p) 

r(2) 



(2) 
iA',, 



i(At) J ro _l(fc|_p) 



fc^_,(Ar)J m+i (^) 



r (2) 

*<7 ~~ 2 



and 



^ fc i;*(-) 

^m.fc T? [d] 



iV (") 6 



(5.27) 



fc_L o-U) 



i(Ar) J m _i{k±p) 



-H 



(i) 
r(l) 



where TV, 



^^^(A^J^x^p) 



is given by Eq. (|5.22| . The Dirac fields 



in the Cartesian and diagonal tetrad representations are 
related via the similarity transformation (|5.9[) . We note 
that S p is given by 1)4. 61| . and S T is defined as 



S T = 



- v /2 



















e-"/ 2 



(5.28) 



We can easily verify that indeed wc have 

1 kj_ ;s 



k v [d] 



(5.29) 
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VI. CONCLUSIONS 

In this paper we have shown step by step how to ob- 
tain the canonical quantization in cylindrical geometry as 
well as in the 1+1 dimensional fluid rapidity and proper 
time coordinates. Assuming that we are interested in a 
quantum back-reaction problem with a background chro- 
moelectric field that is a function only of r, then the sep- 
aration of variables we have made will continue for the 
problem of interest, and the initial value problem choice 
that A(tq) = 0, tells us that the correct mode expan- 
sion at r = To for the Dirac equation is given by the 
problem we have solved. Thus the canonical quantiza- 
tion presented here for the Dirac equation, will allow us 
to determine the initial conditions for the quark state 
(adiabatic vacuum) so that we will be able to numeri- 
cally solve for the production of quarks and the degra- 
dation of the chromoelectric field in a manner similar to 
the problem addressed by Cooper et al. in Ref. [||. By 
additionally having the correct modes for the transverse 
degrees of freedom, we can exploit the cylindrical symme- 
try and reduce the number of momentum modes needed 
in the 3+1 dimensional problems to be TV 2 instead of N 3 . 
These numerical simulations will be soon undertaken in 
the case of color SU(3) quarks. 
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APPENDIX A: COMMENTS ON THE SOLUTION 
OF THE DIRAC EQUATION IN CYLINDRICAL 
COORDINATES 

First, let us note that the same normalization con- 
stant N k ± E is determined by requiring the condi- 
tions [l5j : 

oo 

E ^ti\B\W)*k£mb> 9 ) = *"'> ( A1 ) 

m— — oo 

oo 

E ^rtmb,Q*%&-\B\M = s*- (A2) 

m— — oo 

The following conditions are enforced by construction 

oo 

ni— — oo 
oo 

E ^]- W M*%£ iW M=0- (A4) 



In using the above conditions one makes use of Bcsscl 
functions' property (see discussion regardingthe Graf 
generalization of Neumann's formula in Ref. [22]): 



+ 0O 



m— — oo 



u(x)J m (y)e Hme = e :Bv *J v (z), (A5) 



where 



\/x 2 + 



y 



2xy cos 9 , 



y cos 9 — z cos < 
y sin 9 = z sin < 



(A6) 
(A7) 
(A8) 



Secondly, we discuss the boundary conditions in the 
(x, y) subspace. Provided that we require the solution 
to be regular at the origin, then the generic form (|4.27|) 
reads as follows: for m > 1, we have 



<!>'■ 



A k - e- ie / 2 J m _i(fc ± p) 

2 J m+ i(fc|_/>) 



£l. A k ± e i 



(A9) 



whereas for m < — 1 wc have 



A k± t 
fcj- r<k A 

A_ ^ 

fc-L 4&_L 

A_ n 



iff/2 v 

1 17, 

\6/2 Y 
iff/2 y 

IT. 

Xff/2 y 



m + 



m + 



,(*XP) 



(A10) 



where Y v {x) is the Neumann function. For integer m we 
have the property [22I ] 



Y„ 



J (x) = (-) m J_ m+ ,(x), (All) 
Y m+i (x) = (-r +1 J_ m _ h (x) . (A12) 

Finally, when m = the solution (|4.27[) becomes 



< i;s (p,0) 



A k± e~ w / 2 J_i{k ±P ) 

k±Qk ± e iff/2 Y l(k ± p) 

c k 

A_ 



A k± e i9 / 2 Yi(k^p) 



(A13) 



It is interesting to note that the discussion of the 
boundary conditions in the (x, y) subspace simplifies dra- 
matically if we choose Bessel functions of integer rather 
than half-integer order in (|4.27p . In other words, in the 
(x, y) subspace we can replace the form (|4.27p which in- 
volves Bcsscl functions of half-integer order, and use in- 
stead the form 



^ s (p,9) = e" 



m— — 00 



A k± J m {k ±P ) 
e l6 J m+ i{kj_p) 
C k ± J m (k ±P ) 
A k ^- e w J m+1 (k ±P ) 



A_ ^ 



(A14) 
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which has the correct boundary conditions irrespective of where the form of the solution (|3 . 1 5[) is mandatory [2J 
the value of m. The fact that Eq. (|4.12p allows for a solu- in order to satisfy the orthonormality conditions (|3.19D . 
tion of the form (|A14j) is unlike the case of (t, z) subspace 
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